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Abstract. Under certain constraints on the parameters a, b and c, it is known that 
Schrodinger's equation —d 2 ^/dx 2 + (ax e + bx 4 + cx 2 )ip = Etp, a > with the sextic 
anharmonic oscillator potential is exactly solvable. In this article we show that the 
exact wave function ip is the generating function for a set of orthogonal polynomials 
{Pn\x)} in the energy variable E. Some of the properties of these polynomials arc 
discussed in detail and our analysis reveals scaling and factorization properties that are 
central to quasi-exact solvability. We also prove that this set of orthogonal polynomials 
can be reduced, by means of a simple scaling transformation, to a remarkable class of 
orthogonal polynomials, P n {E) = P„ [E) recently discovered by Bender and Dunne. 
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1. Introduction: 

Recently, Bender and Dunne pQ introduced a remarkable set of orthogonal polynomials 
associated with the one-dimensional Hamiltonian 

H = — f- + x 6 - (4s + 4J - 2)x 2 . (1) 
ax 2 

where J is a positive integer and s = l/4ors = 3/4. The two choices of s correspond 
respectively to the even-parity and odd-parity solutions i/je(%) of the eigen-equation 
Hip = Eip. Bender and Dunne showed that ?Pe{x) is the generating function for a set 
of orthogonal polynomials {P n (E)} in the energy variable E. These polynomials may 
easily be shown to satisfy the three-term recursion relation (with Pq{E) = 1, P\{E) = E) 

P n (E) = EP n ^(E) + lQ(n - l)(n -J- l)(n + 2s- 2)P„_ 2 , n > 2 (2) 

from which it follows that they are orthogonal with respect to a certain weight function 
u(E): 

J P n (E)P k (E)uj(E)dE = 0, n^k (3) 

The 'weight function' u(E) (which we note is not necessarily positive) can be constructed 
by an algebraic method discussed in detail in Ref. |2j. If the initial conditions Po{E) = 1 
and Pi(E) = E, are imposed, each coefficient P n {E) becomes a monic polynomial of 
degree n. The form of the coefficients of the recursion relation satisfied by the polynomial 
system {P n (E)} implies that this system has several remarkable properties. First, the 
squared norms of the polynomials P n {E) vanish for n > J + 1 if J is a positive integer. 
Secondly, each P n (E), with n > J+l, factors into a single product of Jj+i with another 
polynomial, i.e. 

Pj+m+l 

(E) — P J+1 (E)Q m (E), m>0. (4) 

These factor-polynomials Q m (E) form an orthogonal set. There are a number of articles 
devoted to the study the properties of these orthogonal polynomials [2J-IH]- The purpose 
of the present article is to study a set of orthogonal polynomials {P^'(E)} associated 
with the one-dimensional sextic anharmonic oscillator Hamiltonian [TU]-[T%j 

d 2 

H(a, b, c) = — — - + ax 6 + bx 4 + cx 2 , a > 0, (5) 
dx 2 

in which the potential's parameters obey certain constraints. We show that for certain 
constraints on a, b and c, the wave function solution of the Hamiltonian H(a, b, c) is the 
generating function for a set of orthogonal polynomials in the energy variable E. We 
explicity construct the polynomial solvability constraints and prove that they obey a 
three-term recursion relation; consequently, they form a set of orthogonal polynomials 
[TT?] . We study some of the properties of these polynomials such as: for nonnegative 
integer values of J, for which H(a, b, c) is quasi-exactly solvable, the squared norms 
of the polynomials P®(E) vanish for n > J + 1. Further, the polynomials Pff'(E) 
of degree higher than J+l factor into a product of two polynomials, one of which is 
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PjI^E). We also show that under a scale transformation, they leads to the Bender- 
Dunne class of orthogonal polynomials P n (E) = Pf®(E). To this end the paper 
is organized as follows. Section 2 contains a general technique for generating the 
polynomial solvability constraints of the sextic anharmonic oscillator Hamiltonian ([5]); 
these can be easily extended to study the exact solutions for Hamiltonians with even- 
degree polynomials. We further show, through an explicit construction, that the wave 
function solution is a generating function of these polynomials. We prove thereafter 
that this set of polynomials satisfy a three-term recursion relation, and consequently 
they form a class of orthogonal polynomials. In section 4, we show under simple scaling 
transformation the correspondence between a quasi-exact solvable model and a set of 
orthogonal polynomials, and we show that, under suitable change of variables, they 
generalize the class of Bender-Dunne orthogonal polynomials pQ. In the appendix we 
show that some of these polynomials can be expressed in terms of Meixner polynomials 
of the second kind. While the results in the present work including constraining relations 
for the potential parameters, they can usefully be compared with the PT-symmetric 
version of complex sextic potentials, recently studied by Bender and Monou j20j, and 
the work of Bender and Turbiner 



2. Solvability constraints of the sextic anharmonic oscillator Hamiltonian 

Let us assume that the exact solution of the Schrodinger equation 

- ip"(x) + V(x)ip = Eip (6) 

takes the form ip(x) = x(x)e~^ x \ On direct substitution in (JHJ), we obtain the following 
equation for x( x ) 

X "(x) = 2f'(x) X '(x) + (f"(x) - A*) + V(x) - E) X (x). (7) 

Without loss of generality, we may assume, for the node- less eigenstate, that x( x ) is a 
constant, x{ x ) = 1- m this case, Eq.0 reads 

u'{x) = E -V(x) + u 2 (x), (u(x) = f(x)), (8) 

which is a special form of a Riccati equation. For the sextic anharmonic oscillator 
potential V(x) = ax 6 + bx A + ex 2 , we can solve this differential equation exactly for 
certain constraints on the parameters a, b and c. The solvability of this differential 
equation is based on an elegant approach introduced earlier by Rainville [22] providing 
necessary conditions for polynomial solutions of certain Riccati equations. 

Definition 1: By the symbol [yP(x)] ; where P(x) is a polynomial of even-degree, we 

shall mean the polynomial part of the expansion of J ' P(x) in a series of descending 
integral powers ofx. For example: 



Vx 6 - 4x 4 + 7x 2 - 2 



x 3 - 2x. (9) 



With this notation we may state (for a proof, see [22]) 
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Theorem 1: If in 

^ = M*)+u 2 (10) 
Aq(x) is a polynomial of even- degree, then no polynomial other than 

u = ±[^Ao] (11) 

can be a solution of MU\) . If the degree of Aq is odd, there is no polynomial solution 
°f 11U\) - As an example, for the first-order nonlinear differential equation v! = 
2 — 7x 2 + Ax 4 — x 6 + u 2 , we have u(x) = — (x 3 — 2x) is a solution which fact can 
be easily verify through direct substitution. 

By means of this theorem, we can search for exact solutions of the differential equation 
(JHJ) with V(x) = ax 6 + bx 4 + cx 2 , in which case we have 

±(v^x 3 + -^x), (12) 



u = ± 



V ax & + bx 4 + cx 2 — E 



if 12a 2 — b 2 + 4ac = and E = For a physical acceptable solution satisfying 

ip(±oo) = 0, we have, for u(x) = f'(x), that 



f(x) = ^x 4 + -^x 2 (13) 
4 Ay/a 

Consequently, for the Schrodinger equation 

- ip"(x) + (ax 6 + bx 4 + cx 2 )ip = Eip, (14) 
we may assume that exact solution takes the form 

ip(x) = x( x ) e 4 4VE (15) 

which has been adopted in the literature for this class of potentials PUJ-JEI- in order 
to find x( x )> we notice by means of (|15p and (|7J) that 

X " = 2(v^x 3 + -^x) X ' + -E+ (3y/E - ^ + c)x 2 ) X (16) 

which clearly yields a ground-state eigenenergy = if the potential parameters 
satisfy the constraint 3a/o" — j- + c = 0. The search for polynomial solutions x( x ) = 
S of (fTfij) can then be made by means of the standard techniques of series solution 
of second-order differential equations. However, these exact solutions can be explicitly 
generated by means of the asymptotic iteration method (AIM) recently introduced [23J. 
Actually, we should stress the usefulness of AIM as a method for determining the explicit 
form of the solvability constraint polynomial P n (E). AIM was first introduced |2Sj to 
solve second-order linear differential equation of the form 

y" = X (x)y' + s (x)y (17) 

where Xq(x) ^ and sq(x) are sufficiently many times continuously differentiable. 
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Theorem 2: Given that Xo(x) 7^ and Sq(x) are sufficiently many times continuously 
differentiable, the second-order differential equation has the general solution 



u\x) 



exp (- J adt) C 2 + C x J exp ( J (A (r) + 2a(T))drjdt 



(18) 



a (19) 



if for some n > 0, 

"Sn "Sn— 1 
An A n _i 

where 

— Xn-l + S n-1 + A()A n _l, and S n = s' n _ 1 + SoA n _i- 

The asymptotic iteration method was soon adopted [2II|-[2I3 to investigate the solutions 
of eigenvalue problems of Schrodinger type. In such applications, one immediately faces 
the problem of transforming the Schrodinger equation (with no first derivative) into 
the form ([170. The use of asymptotic solutions of the Schrodinger equation under 
consideration is the usual approach employed to overcome this problem. It is important 
to mention that the asymptotic form is very crucial for the convergence of the iteration 
method to exact solutions. For the sextic anharmonic oscillator potential V(x), or for 
more general even-degree polynomials V(x) = Ya^i 1 a i x2 \ n — 2, the construction 
based on Rainville's approach [22] provides a straightforward technique to generate a 
proper asymptotic form that stabilizes and accelerates the convergence of AIM. The first 
few iterations of (fTTjj) yields polynomial expressions {P n (E)}. The exact eigenvalues can 
be computed, in turn, as the zeros of these polynomials (for convenience, we denote 
P (E) = 1,P 1 (E) = 1): 

• P 2 (E) = 1Ja\E - b, if 12a 3 / 2 - b 2 + Aac = 0, x 2 (z) = 1. 

• P 3 (E) = 2^/EE - 36, if 20a 3 / 2 - b 2 + Aac = 0, xs{x) = x. 

• P 4 (E) = AaE 2 - Yl^fabE + 24a 3 / 2 + 36 2 + 8ac, if 28a 3 / 2 - b 2 + 4ac = 0, 

Xt(x) = P (E) - ^-x 2 . 

Ay/ a 

• P 5 (E) = AaE 2 - 20^/a~bE + 120a 3 / 2 + 156 2 + 24ac, if 36a 3 / 2 - b 2 + Aac = 0, 

X5 (x)=x(P 1 (E)-^lx 2 ). 



P Q {E) = 8a 3 / 2 E 3 -60a6E 2 + (720a 2 -90 v /^6 2 + 112a 3 / 2 c) J B-552a 3 / 2 6-156 3 -120a6c, 
if 44a 3/2 - b 2 + Aac = 0, 

Xe{x) - P (E) - —x 

P 7 (E) = So?/ 2 E 3 - 8AabE 2 + (1680a 2 + 2Wy/ab 2 + 208a 3 / 2 c)£ - 3480a 3 / 2 6 - 1056 3 - 
5QAabc, if 52a 3 / 2 - b 2 + Aac = 0, 
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It is quite clear that the even-parity wave function solutions of sextic oscillator 
Hamiltonian (jSJ) satisfy 

^«W = g^gf^. -0,1,2,..., (20, 
while for the odd-parity wave function solutions we have 

By means of the differential equation (jUJ), we see that the polynomial solvability 
constraints satisfy the recursion relations: 

• For even-parity solution, with Po{E) = 1, 

P 2n .+ 2 (P) = {2y/EE - {An + l)b)P 2n {E) 

+ 2n{2n - l)[4(4n - 1)J - b 2 + Aac]P 2n _ 2 {E) , (22) 

for n = 0, 1, 2, ... , (Note that for n = the nP 2n _ 2 {E) term is not present) subject 
to the condition 

A{An + 3)a§ - 6 2 + 4ac = (23) 

• For odd-parity solution, with V\{E) = 1, 

V 2n+ 3(E) = {2^E - {An + 3)b)V 2n+1 {E) 

+ 2n(2n + l)[4(4n + l)at - 6 2 + 4ac]P 2n _i(P), (24) 
for n = 0, 1, 2, ... , (for n = the nV 2n -i{E) term is not present) subject to 

A{An + 5)a§ - 6 2 + 4ac = 0. (25) 

3. Quasi-exact solvable systems and orthogonal polynomials 

The Hamiltonian (jSJ) has the following scale transformation property 

H{a, b, c) = a*H(l, ba~%, ca~%) (26) 

Setting ba~i = 2t and ca~^ = t 2 — AJ — 3, where J is a nonnegative integer and t is a 
real number, the Hamiltonian (jSJ) reads 

P = — + x 6 + 2tx A + {t 2 - A J - 3)x 2 . (27) 
ax 2 

In this case, P n {E),n = 0, 1, 2, . . . becomes 

P 2 n(£(a,M) = P ^E {iba _,^ h) ) = P n {a^E) = 2"a"P«(P), n = 0, 1 

P 2n -2{E (aM ) = P 2n _ 2 {a^E {iba _ lca _ h) ) = P n - X {a*E) = Z^a^ P^E), n = 1, 2 
^ P2n+2{E {aM} ) = P 2n+2 {a^E {i ba _ i ca _ h) ) = P n+1 {a^E) = 2 n+1 a^ P^ +1 {E) n = 0, 1 

and the recursion relation ()22|) now reads for n = 1,2,... 

P^P) = (J5 - (4n + l)t)pW(^) + 8n(2n — l)(n — J — l)Pil(P), (28) 
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which uniquely determines, with P^{E) = 1 and P?{E) = (E - 5t)P (t) (P), all the 
polynomials P®(E),n = 1,2, ... in terms of P (P). With these initial condition, the 
recursion relation ()28|) generate a set of polynomials, the next four of which are 

p[ t) {E) = E-t 

P 2 W (P) = E 2 - QtE + 5t 2 - 8 J 

P 3 (t) (P) = P 3 - IhtE 2 + (48(1 - J) - 8 J + 59t 2 )P - 3t(16 - 40 J + 15t 2 ) 
P 4 (t) (P) = P 4 - 28tP 3 + (288 - 176* + 254t 2 )P 2 - 4t(528 - 392 J + 203t 2 )E 
- 48(45J-38)t 2 + 585t 4 

On multiplying the recursion relation (|28)l by E n u{E) and integrating with respect to 
E, using the fact that P^\E) is orthogonal to P fc , k < n, we obtain a simple, two-term 
recursion relation for the squared norm j n of Pjp as: 

7n = 8n(2n-l)(J-n + l)7 n _i, (29) 

which is independent of t. The solution to this equation with 70 = 1 is 

7n = ft 8*(2* - 1)(J - k + 1) = ^(T+J-r!)^ - (30) 

The interesting factorization property pQ follows when J takes nonnegative integer 
values. This is clear because the third term in the recursion relation (f2*Kj) vanishes 
when n = J + 1, so that all subsequent polynomials have the common factor Pj + i(E). 
To illustrate this factorization, we list in factored form the first five polynomials for the 
case J = 1 

PV(E) = 1, 
P?{E) = E-t, 
P?{E) = E 2 - QtE + 5t 2 - 8 
P 3 W (P) = (E- 9t)(E 2 - QtE + 5t 2 - 8) 
P 4 W (P) = (P 2 - 22tP + 117t 2 + 120) (P 2 - QtE + 5t 2 - 8) 
P 5 (t) (P) = (P 3 - 39tP 2 + (568 + 491t 2 )P - 1989t 3 - 6072t)(P 2 - QtE + 5t 2 - 8) 
In general 

P^jAE) = Q^{E)pfl{E) (31) 

Substituting (J3T| into (|28jl . one can obtain the recurrence relation immediately for the 
factor-polynomial Q$(E): 

QW(P) = (P-(4n+4J+l)t)gi*! 1 (P)+8(n-l)(n+J)(2n+2J-l)Qjl 2 (P)(32) 

with initial condition Qq\e) = 1 so that Q$(E),n = 0, 1,2, . . . are again orthogonal 
polynomials. The squared norm of Q$(E) is given by (with 7^ = 1) 

n 

In = II 8k ( k + J+l)(2k + 2J+ 1). 
k=l 
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Further, for the odd parity case, the recursion relation (|24|) reads, under the scaling 

dm, 

VZUE) = {E- (An + 3)t)P®(E) + 8n(2n + l)(n - J - \)V { ^E). (33) 
with Vq\e) = 1. The first few explicit polynomials are 
V?(E) = E-3t 

P£\E) = E 2 - lOtE + 21t 2 - 24 J 

7?W (£) = E 3 — 21tE 2 + (80 - 104 J + 131t 2 )£ - 3t (80 - 168 J + 77t 2 ) 

Pi t} (E) = £ 4 - 36tE 3 - 12tE(400 - 312 J + 183t 2 ) + £ 2 (416 - 272 J + 446t 2 ) 

+ 9(448(J - 2) J - 16(-74 + 77 J> 2 + 385t 4 ) 

Again on multiplying (|31|) by E n ~ 1 oj(E) and integrating with respect to £7, using the 
fact that V^'(E) is orthogonal to E k , k < n, we obtain a simple, two-term recursion 
relation for the squared norm 7^: 

^ = 8n(2n + l)(n- J - l) 7 „-i, (34) 

The solution to this equation with 7 = 1 is 



n 



k=i L \ J n + X J 

It is clear that the squared norms (f37I|) and (p$H|) vanish for n > J + 1, as expected. The 
classes of orthogonal polynomials discovered by Bender and Dunne follow directly by 
setting t = in (j2H| and (jHHJ). The factorization property for the polynomials {V!jp(E)} 
in the case of J = 1 can be illustrated by means of the polynomials 

V$\E) = 1 
V { i\E) = E-3t 
V?(E) — (E — 7t)(E - St) - 24 
ri t] (E) = ((E- 7t)(E - 3t) - 24) (E - lit) 
V^(E) = (168 + (E — lbt)(E - llt))((E - 7t)(E - 3t) - 24) 
Vi t] {E) = ((E - 7t)(E - 3t) - 24) (E 3 - AhtE 2 + (744 + 659t 2 )£ - 3135t 3 - 9528t) 
In general 

rHl J+1 (E) = Q^\E)vf +1 {E) 

where Q n (E) satisfy 

Qt ] {E) = (E-(4n + 4J + 3)t)Q { S 1 (E) + 8(n-l)(n + J)(2n + 2J+l)Q<il 2 (E), n>l. 
The squared norm for these polynomials is then 

n 

In = II 8k (k + J + l)(2fc + 2J + 3) 
fc=i 

Clearly, the squared norms of the polynomials Q$ (E) and Q$ (E) do not vanish. The 
weight functions for all of (E) and P« (E) as well as of the polynomials Q%> (E) and 
Q$(E) can be computed by means of the method discussed in 
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4. Conclusion 

In the field of quantum mechanics the antecedent to the concept of quasi exact solutions 
may perhaps be found in an early paper by Wigner in 1929 in which the following 
simple idea is explored: first choose a wave function, and then find the corresponding 
potential. The form originally chosen for the wave function was the exponential of a 
polynomial, and it was shown that the construction worked out if the coefficients met 
certain conditions. As is evidenced by the references in this paper, and indeed in the 
contents of the paper, many results have been discovered since that early work of Wigner. 
In order to test general approximation theories (so essential to applications of quantum 
mechanics), it is always useful to have at hand a large collection of exact solutions. 
This is one clear area of utility for the outcome of this work. Of course, problems that 
start as questions in physics soon take on a life of their own and may later generate a 
repository of results that serve the field from which they originally emerged. 
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Appendix: 

In this appendix, we search for closed-form expressions for the orthogonal polynomials, 
{P$} defined by (j2*8Jl and {V®} defined by (JHHJ), in terms of hypergeometric orthogonal 
polynomials |2Zj- We consider first the case of polynomials defined by the recurrence 
relation (|28|). Since its squared norm (|30p vanishes for n > 1 + J, we can take 
J = n + i, 2 = 0,1,2,.... Thus we write (J25J) as 

P^UE) = {E- {An + l)t)P£\E) - 8n(2n -l)(i + (36) 

This recurrence formula can be compared with 

P n+1 (E) = (E- (dn + f))P n (E) - n(gn + h)P n ^(E). (37) 

studied in [2E] With a 2 = Ag - d 2 > 0, 5 = ±, r] = 1 + | and the choice 2f = Srja, the 
recurrence formula (|37|) becomes, for n > 0, 

M n+1 (E) = [E — (2n + n)5]M n {E) - (5 2 + l)n(n + n - l)M n -\(E) (38) 

where 

M n (E;S,r ] )= Q"p„^) (39) 

is a Meixner polynomial of the second kind. For our case, we have d — At, f — t, g — 
16 (i + 1), h — —8(i + 1), and, using the notation 



a 2 =Ag-d 2 = 16(Ai + A-t 2 ), or a = A^Ai + A - t 2 
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where 



(2VM + 4 - t 2 E^j . (42) 



and 

. d t , h 1 

£ = - = — ; n = l + - = - 

the recurrence formula (j3ti)) becomes 

M n+1 (£) = [E [ 4w + 1}t ]M n (E) - 2n{2n + —M n _i(E), (40) 

+ V ; 1 2V4z + 4-t 2J V ; V4i + 4 - t 2 V ; V ; 

MJE) = MJE; ; t -) (41) 

1 

,2V4i + 4-t 2 , 

For the orthogonal polynomials {V%'(E)} defined by (p^jl . we again let J = n + i, i = 
0, 1,2,.. .. Thus we write (jSHj) as 

= (£7 - (4n + 3)t)V®{E) - 8n(2n + l)(i + \)V® X (E), (43) 

which can be compared, again, with (j!T7jl for d — 4t, f — 3t, g — 16(i + 1), /i = 8(i + 1). 
In this case, we have 

a 2 = Ag-d 2 = 16(4i + 4 - t 2 ), or <r = 4V4i + 4 - t 2 

and 

d t h 3 

(5 = - = — r7 = l + - = -, 

<r V4i + 4-t 2 ' 5 2' 

and the recurrence formula ()43j) becomes 

M„ +1 (£) = [E — ^ 4n + 3)t ]M^)- 2 ^ + 1) ^ + 1) M w _ l( E),(44) 
+U ; L 2V42 + 4-t 2J 1 ; V4i + 4 - 1 2 U j ' 1 J 

MJE) = MJE; . = -) (45) 



where 



1 



2V4i + 4 - t 2 



^2V4z + 4 - t 2 £^ . (46) 
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